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Preface 


The enclosed report is a copy of a paper submitted for publication 
in the Journal of Composite Materials . The Issuing of this report is 
not intended to duplicate the Journal’s publication effort. This report 
is intended to accelerate the dissemenation of information resulting 
from investigations of advanced composite materials. 


Abstract 


The cured shape of unsymmetric laminates do not always conform to 
the predictions of classical lamination theory. Classical lamination 
theory predicts the room-temperature shapes of all unsymmetric laminates 
to be a saddle. Experimental observations, however, indicate some 
unsymmetric laminates have cylindrical room-temperature shapes. In 
addition, some unsymmetric laminates, exhibit two stable room-temp- 
erature configurations, both cylindrical. This paper presents a theory 
which explains these charo/jteristlcs. The theory is based on an exten- 
sion of classical lamlnacion theory which accounts for geometric nonlin- 
earities. A Rayleigh-'Aitz approach to minimizing the total potential 
energy is used to obcain quantitative information regarding the room- 
temperature shapes of square T300/5208 [0,2/902^1' B^^aphite- 

epoxy laminates. It is shown that, depending on the thickness of the 
laminate and eho length of the side of the square, the saddle shape 
configuration is actually unstable. For values of length and thickness 
that render the saddle shape unstable, it is shown that two stable 
cylindrical shapes exist. The predictions of the theory are compared 
with existing experimental data. 


Introduction 


Most calculations used to predict the response of laminates to 
static, dynamic, or thermal loadings are based on what has become to be 
known as classical lamination theory [1], [2], [3]. Tills is a linear 
theory and Is based on the following major assumptions; 

(1) the displacements are continuous throughout the laminate, 

(2) the Kirchhoff hypothesis regarding undeformed normals is 
assumed to be valid, 

(3) the strain-displacement relationship is linear, 

(4) the material is linearly elastic, and 

(5) the through-the-thlckness stresses are small in comparison to 
the in-plane stresses. 

The theory smears the individual lamina properties by integrating 
the constitu'''.ive equations through the thickness of the laminate. As a 
result of this integration, force and moment resultants are defined. In 
addition, the well-known A, B, and D matrlcles are defined. Wliile this 
theory is quite capable of accurately predicting static deflections, 
natural vibration frequencies and mode shapes, buckling loads and mode 
shapes, and thermal expansion coefficients of laminates, there are 
physical situations for which the theory fails to predict the correct 
answer. Two situations of note are: the Inability of the theory to 

predict the response of thicker laminates, and; its inability to explain 
the behavior of laminates near edges. The former problem has been 
studied by several investigators [4], [5} and satisfactory corrections 
to the theory have been obtained. The latter problem Is now a classic 
and has been studied by many individuals, at least for the case of the 
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straight free-edge, K survey ef the edge problem has been put forth in 

[ 6 ]. 

There appears to be another situation for which classical lami- 
nation theory fails to give the correct answer. Specifically, it 
appears the theory is unable to correctly predict the room- temperature 
shapes of thin unsymmotric laminates. Hyer [7] has documented the room- 
temperature shapes of several families of unsymmetric laminates and 
found that the room-temperature shapes of some thin unsymmetric lami- 
nates are closely approximated by right circular cylinders. In addi- 
tion, some thin laminates have two room-temperature cylindrical shapes. 
These results are in contrast to the predictions of the classical the- 
ory. The classical theory predicts the room-temperature shapes of all 
unsymmetric laminates to be- a saddle with unique (slngle^valued) curva^ 
tore characteristics. Specifically, Hyer found that 100 x 100 mm and 
150 X 150 mm [Og/QO^l.^, T300/5208 graphite-epoxy laminates cured to 
become cylindrical at room temperature. In addition, they exhibited a 
snap-through or oil-canning phenomenon. The cylindrical shape could be 
snapped into another cylindrical shape which had the same characteris- 
tics as the first shape. However, the second cylinder was oriented 
perpendicular to the first cylinder and the curvature of the second 
cylinder was of opposite sign. Hyer showed that thicker (say, 8-layer) 
100 X 100 nmi unsymmetric laminates conformed to the predictions of the 
theory. 

To explain the behavior of unsymmetric laminates, it was assumed 
that it would be necessary to incorporate a nonlinear effect into clas- 
sical lamination theory. The existence of two room-temperature shapes 
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(i.e. the two cylinders) essentially ruled out a linear extension to the 
theory since a linear extension would lead to the prediction of a 
unique shape, albiet perhaps not a saddle shape, Fur thermo re, since the 
out-of-plane deflections pf the u.tsymmetric laminates were on the order 
of many laminate thicknesses, geometric nonlinearities ware felt to be 

ft " 

an important effect, Tims, in an effort to explain the behavior of 
unsymmetric lafirlnates, classical lamination theory was extended to 
Include geometric nonlinearities through the strain-displacement rela- 
tionship. This extension was applied to the analysis of the tOj^/90^)^, 
n“l,2,,.. family of laminates. This family was chosen for study because 
this class of laminates always exhibits a snap-through phenomenon, there 
is data available for different thickness (0 /90 laminates, and due 
to some, of the , B^.j, and terras being zero with this family, the 
algebra associated with this family of laminates is simpller than the 
algebra associated with other families. The analysis of the room- 
temperature shapes of this family of laminates is the subject of this 
paper. The paper traces the development of the analysis which success- 
fully predicts the existence of two room-temperature cylindrical shapes, 
and compares the predictions with the available data. 


Problem Formulation 

Since a problem formulation which includes geometric nonlinearities 
would result in nonlinear governing equations, it was assumed from the 
beginning that obtaining a closed-form exact solution for the unsym- 
raetric laminate problem would be difficult and not really necessary. 

The occurence of the cylindrical shape is so prevalant with thin lami- 
nates that it was hypothesised that one is dealing with a fundamental 
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phQnomanon rather than sem® higher order effect. Tlius any good approx- 
imate theory would reveal the mechanica of the problem. Tlie prol4em is 
idealized as follows. A cur<;d laminate (and uncured prepreg) is flat at 
the elevated curing temperature, fig. la. As the laminate cools, it is 
assumed it is free from any external mechanical forces which produce net 
work. It is assumed the out-of-plane deflections develop only because 
of the differences in thermal expansion properties of the Individual 
lamina, Tills Idealization ignores the effects of any mechanical con- 
straints the autoclaving and vacuum bagging process may exert on the 
laminate. Upon cooling, the laminate deforms into one of the shapes 
given by figs, lb, c and d. Of these shapes, the one that actually 
occurs Is the one associated with a minimum of the total potential 
energy. The shape scenario given by figs, lb, c and d includes the 
saddle shape observed for the thicker laminates and predicted by the 
clt:.iiical theory, and the two possible cylindrical shapes observed for 
the thinner laminates. It is assumed that in attaining these shapes, 
geometric nonlinearities are important. 

Since it is assumed external tractions are not important during the 
cooling process, the total potential energy, including the effects of 


thermal expansion, is given by [8] , 


W 


WdVol, 


Vol 


^ “ 2 ^ijkl®ij®kl 




( 1 ) 

( 2 ) 


where W equals the strain energy density. The are the elastic 

constants of the material and the 8^^ are coefficients related to the 
elastic constants and the coefficients of thermal expansion of the mat- 
erial. Both the elastic properties and the thermal expansion coeffi- 
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cients are aaauraed to be tenipetature-independent, The are the 
strains In the material and AT Is the temperature change in the material 
due to cooling from curing. In eq. (2), since the problem is a plane- 
stress formulation, i and J assume the values 1 and 2. These values are 
not directly related to the principal material directions of the lamina, 
but rather, 1 and 2 are associated with the x and y directions of the 
laminates (see fig. 1). Thus the following relations apply; 
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( 6 ) 

(7) 

( 8 ) 


As usual, z»0 is the midplane of the laminate. The quantity u° is 
the laminate midplane displacement in the x-direction, v® is the larai- 
nafe midplane displacement in the y-direction, and w is the out-of-plane 
displacement of the midplane. 

Equations (6) - (8) represent the principal departure from clas- 
sical lamination theory and include the usual approximations associated 
with thin-plate theory when employing the nonlinear geometric effects in 
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the stmln-displacdment: relations. These approxlmatlona assume the 
elongation and shearing strains and the squares of the rotations are the 
same order of magnitude and this order is small compared to unity [9], 
[103. 


For a laminate the can be related to the reduced 

stlffuesses, and the can he related to Che Q^j's and ot^, a^, and 

Ot^y, the laminate thermal expansion coefficients in the x-y coordinate 
system. Expanding eq. 2 yields 


W 


i ^11®11 '*■ %2® ll ‘*22 ^^11®12 I ^22*^22 


“ ^%l“x ^12“y^®ll^’^’ " ^^12“x ^22^y^®22^'^ 


(9) 


(a “O for this family of laminates.) Tlie problem has now been reduced 
xy 

to one of finding the deformation u®, v“, and w as functions of x and y 
which, through eqs. (2) - (9), minimize eq. (1). As previously men- 
tioned, approximate solutions Co u“, v®, and w are sought. In seeking 
realistic approximate solutions, two basic assumptions are made. First, 
It is assumed that even in attaining the cylindrical shape, the midplane 

elongation strains, e® and e®, do not vary much from the linear predic- 

X y 

tion (l.e. e® and e" independent of x and y) . Second, it is assumed 
X y 

that to the order of the nonlinearity considered here, the midplane 
shear strains are negligible i.e, e^® « 0, Since for lami- 

nates classical lamination theory predicts e ® to be zero and e® and 

xy X 

e° to be constant, these two assumptions could be lumped into one by 
saying it is assumed that even in attaining the cylindrical shapes, 
magnitude of the midplane strains do not vary much from the predictions 
of the classical theory. However, for rationalizing the choice of the 
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functional form of tha approximate aolutlons, the two laBuea are eon- 
arated* 

It is asaumed w(x,y) ie of the form 

w(x,y) - Y (nx^ + by^) , ( 10 ) 


a and b being conatanta. With this functional form for w, both the 

claBsical lamination solution, a * -b, and either of the two cylindrical 

shapes, fig, Ic and d, can be approximated* For fig. Ic the solution 

can be a 0, b » 0 while for fig. Id the solution can be a * 0, b 0. 

Using the kinematic assumptions regarding the midplane strains, e®, 

e°, and e ®, the approxlmaie solutions for u® and v° are given by 
y 


2 3 , 2 

) * cx - - MEL. 


( 11 ) 

( 12 ) 


e and d being constants. Using eqs. (H) and (12) in eqns. (6) - (8) 
yields 



(13) 

, abx^ 
- d - ^ 

(14) 

"0. 

(15) 


Note that if It were not required to have e^“ be zero, the third term in 
each of eqs. (11) and (12) would not need to be included and then the 
second term in eqs. (13) and (14) would not appear. However, with the 
[0^/90^]^ family, it is felt In-plane shear strains are Impossible and 
this is the prime factor in choosing the functional form for u° and v®. 

At this point a, b, c and d are considered as generalized coordi- 
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nates and are to be determined. The problem of finding a minimum for 
the total potential energy, W, becomes a problem of finding solutions to 
the values of a, b, c, and d so that the first variation of W is zero, 
i.e. 

6 W "> (|“)da 4* (|^) 6 b 4* 4- (|j)^d 5 0 . (16) 

Tills variation is done with all the assumptions in the Introduction 
with the exception of (3), 

Calculations Associated with the Solution 

The X'-y-z coordinate system in fig. 1 is assumed to be situated 
such that at the elevated curing temperature the laminate is defined by 
the region 

-L /2 < X < L /2 

A A 

-L„/2 < y < L /2 (17) 

tX / 

-h /2 < z < h /2 . 

With these limits on the spatial variables and with the various material 
properties involved, eq. ( 1 ) takes the form 
L /2 L /2 h /2 

W- W(a,b,c,d,Q^j,a^,ay.AT,x,y,z)dxdydz. (18) 

X--L /2 y»-L /2 z“-h/2 
A y 

The Involved, but straight-forward, process of substituting eqs. (10) 
and (13) - (15) into eqs. (3) - (5), substituting these results into eq. 
( 9 ), performing the spatial integrations in eq. (18), and finally taking 
the first variation, eq. (16), leads to an equation of the form; 

6 W « fj^(a,b,c,d) 6 a 4- f 2 (a,b,c,d) 6 b 4- 


f 2 (a,b,c,d) 6 c 4- f^(a,b,c,d)(Sd » 0 . 


( 19 ) 
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Equaclon (19) Immediately leade to four equatlonas 

2 

f j^(a I b| C)d) • ~ Cj^cb Hf* C2ab + ZC^ab 4* 


D^^a - C^cb + ZCgab^ - Cgdb + 

Dj^2b « C^db + Cgab^ + Cgb^ + 

aJ/48)N]Jb + + (L^M8)Njb - 0 , 

2 2 2 

fgCajb.Cjd) - -Cj^ac + C2a b + ZC^a - C^ac + ZC^a b + Dj^2a 


( 20 ) 


Cgds ^ Cydfii 4* Cg£i b Hh 2C|^db ^ B22^^ ** 

Dazb + (by/48)NjJa + (L^/48)Nja + mJ - 0 , 


(21) 


f2(a,b,c,d) ■ A^j^c - Cj^ab - + A^^2*^ *■ ^4®^ “ ^ ® * 


(2Z) 


,T - 0 . 


(23) 


f^(a,b,c,d) - A^ 2 ° " ~ ®2Z^ 

The conatants are defined In the Appendix and B^^ and 

have the familiar definitions associated with laminates, The other 
definitions used in eqs, (20) - (23) are, 
h/2 

/ - AT J (Q^^a^ + ^2«y>‘^'' » ^24) 


-h/2 

h/2 
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- AT 

y 


-h/2 
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( 26 ) 
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M. 



12^x • 


(27) 


These quantities are immediately recognizable as the effective in-plane 

thermal loads, and nL and the effective thermal moments, m 3 and 

X y A y 

It should be noted that when L^»Ly*0, the coefficients through are 
all zero and eqs. (20) - (23) reduce to the equations of classical lami- 
nation theory. 


Solution of Equations. Numerical Results 

Solutions to eqs, (20) - (23) were obtained by solving eqs. (22) and 
(23) for c and d in terms of a and b and substituting these relations 
into eqs. (20) and (21). Thus eqs. (20) and (21) become coupled cubic 
equations for the quantitites a and b. These two resulting equations 
have the characteristics of being able to be reduced to a single cubic 
equation for either a or b. Such an equation would have either one or 
three real roots. This reduction approach was not used, however, and 
eqs. (20) and (21), in terms of a and b, were solved numerically. Solu- 
tions were obtained for several laminates using elastic and thermal 
expansion properties of T300/5208 graphite-epoxy. It was assumed the 
curing temperature of T300/5208 is 177“ C (350° F) and that the lami- 
nates are cooled to a room temperature of 21“ C (70“ F) . The material 
properties used in the calculations were: 

= 181 GPa (26.2 x 10^ psi) 

= 10.3 GPa (1.49 x 10® psi) 

^12 “ 

Gj ^2 ■= 7.17 GPa (1.04 x 10® psi) 
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Ct^ « -0.106 X 10"^/ » C C-0.059 X 10”^/* F) 

«2 25.6 X 10"^/“ C (14,2 X lO"^/* F) 

The elastic properties were taken from (2] while the thermal expansion 
coefficients were taken from [11] , Solutions were obtained for square 
laminates (L «L «L) ranging from 0 to 150 nm li,v length on a side. Two 

* y 

thicknesses were considered, £02^^®2^T data 

were available for some checking of these two thickness cases. Figure 2 
shows the characteristics of the predicted room-temperature shapes of 
the t02^^*^2^T 3 shows the characteristics for the 

thicker laminate, Each figure Illustrates the effect of the size of the 
laminate, L, on the room temperature shapes. 

Immediately obvious from the figures is the existence of t hree 
possible room-temperature shapes of the laminate if the lengths of the 
sides are greater than some critical value. For both l^ii*, „nates , at L«0 
the room-temperature shape is the saddle predicted by classical lamina- 
tion theory, a^-b. This solution is denoted by point A on figs. 2 and 
3. As the sides of the laminate increase in length, say to L«=25 mm for 
the [ 02 / 902 ],!, laminate of fig. 2, the shape is still predicted to be a 
saddle but one which is shallower than the one predicted by the clas- 
sical theory. As the lengths of the sides increase, the saddle shape is 
still predicted to exist but it gets shallower and shallower. At some 
critical length, the solution bifurcates. For the thinner laminate the 
critical length is 35 mm while for the thicker laminate the critical 
length is 71 mm. The bifurcation point is denoted as B on the figures. 
For lengths greater than the critical length, three room temperature 
shapes, each represented by a different branch on the figures, can 



12 


possibly exist. These branches are denoted as BC, BD, and BE on the 
figures, Branch BD represents a continuation of the saddle shape 
(a^-b) but the other two branches represent a radical departure from a 
saddle shape. Branch BC represents a shape which has a large curvature 
In the x-directlon and practically no curvature In the y-direction, fig. 
Ic. On the other hand, branch BE represents a shape which has a large 
curvature In the y-direction and very little curvature in the x-direc- 
tlon, fig. Id. The shapes associated with these latter branches can be 
considered cylindrical because as the laminate gets larger, i.e. L 
increases, the one curvature asymptotically approaches zero while the 
other curvature asymptotically approaches a non-zero constant value. As 
seen from the figures, the latter two branches have certain symmetry 
characteristics. These symmetry characteristics are such that for a 
given length, the values of a and b associated with branch BC are equal, 
respectively, to the values of -b and -a associated with branch BE. 

Figures 2 and 3 show both laminates exiblt similar behavior. There 
are differences however. The two main differences are that the curva- 
tures for the thicker laminate are less than the curvatures for the 
thinner laminate, and, the critical length for the thicker laminate is 
greater. Thus, compared to a [ 02 / 902 ]^ laminate, a [0^/90^]^ can be 
made larger before the triple-shape phenomenon occurs. 

Stability of the Predicted Shapes 

With multiple solutions to a nonlinear problem, the question arises 
as to the stability of the various solutions. If any of the solutions 
do not represent a stable solution, those solutions will not be physi- 
cally realizable. Equating the first variation of the total potential 
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energy to zero yields equilibrium positions of the laminate which either 
maximize or minimize the total energy. For stable equilibrium, the 
total potential energy must be mimimlzed. Thus, for stable equilibrium 
the second variation of the total potential energy, 6 W, must be pos- 
itive definite. Fvom stability theory [12] for this discretized system, 
stability of the tbrium positions for the laminate is possible if 

and only if the following matrix of coefficients is positive definite; 

9a 9b 9c 9d 

9f2 9 f 2 ^^2 ^^2 

9a 9b 9c 9d 

S 

9a 9b 9c 9d 

9f, 9f, 9f, 9f, 

4 ^ ^ 4 

9a 9b 9c 9d 

For a given equilibrium solution at a given length, e.g. the saddle 
solution of the triple-valued solution at L=100 mm, each element of the 
matrix is evaluated numerically by substituting in the values of a, b, 
c, and d corresponding to that solution. If each of the principal 
minors of the matrix are positive definite, the matrix is positive 
definite and the solution corresponds to a stable equilibrium solution. 
Otherwise the solution corresponds to an unstable equilibrium solution. 
Using this scheme for the solutions shown in figs. 2 and 3, it was found 
that the saddle solutions corresponding to the single-valued solutions, 
segments AB in the figures, were stable . On the other hand, the saddle 
solutions corresponding to the triple-valued solutions, segments BC in 
the figures, were unstable . The other two branches of the triple- 
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valued solutions, segments BD and BE, represent stable solutions. 
Physically this means that for square laminates of the 
family, if the length of the sides of the laminate exceed 35 mm, the 
saddle shape does not exist. Instead, two cylindrically shaped equili- 
brium configurations exist. If the laminate is from the thicker 
[0^/90^^]^ family, the length of a side must exceed 71 mm before the 
saddle-shape equilibruim configuration disappears and dual cylindrical 
shapes appear. The fact that two stable cylindrical equilibrium solu- 
tions are predicted to exist is felt to be significant since it corre- 
lates well with the the reported snap-through phenomena associated with 
these types of laminates. 

Experimental Results 

Shown on fig. 2 are two data points. These points correspond to 
the curvatures of cylindrical measured by Hyer. 

Figure 3 shows one data point. This point corresponds to the curvature 
of a [+45^/-45^]^ saddle-shaped T300/5208 graphite-epoxy laminate as 
measured by Pagano and Hahn [13] . The comparison between magnitudes of 
the predicted and experimentally measured curvatures is fair. More 
Importantly, however, the character of the measured shapes, i.e. cylin- 
drical or saddle, compares well with Che predictions. For the 150 x 150 
mm laminate shown in fig. 2, only the major curvature of 11 laminates 
were measured in the original work and the average curvature and the 
range are shown here arbitrarily as a y-directlon curvature, b. These 
curvatures could have just as easily been called an x-direction curva- 
ture. In this case the experimental data would have been associated 
with the variable a. Also the major curvature from the one 100 x 100 mm 
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specimen Is also arbitrarily associated with the y-dlrectlon curvature. 
For the data point in fig. 3, the curvature of the specimen was never 
measured directly. The out~of-plane deflection across the diagonals of 
a 63.5 mm (2.50 in.) square laminate was measured and the curvature was 
computed from this measure. Again this single curvature measurement was 
arbitrarily associated with the variable b. 

It should be noted in fig. 2 that the x-direction curvature for the 
100 X 100 mm laminate was measured to be slightly negative. The theory 
predicts this curvature to be slightly positive. The reason for the 
descrepancy is not clear. However, it is not felt to be due to mea- 
suring error in the experimental determination of curvature. This 
points needs further investigation. 

Discussion 

Despite the lack of large amounts of quantitative experimental data 
to compare with the complete range of the numerical predictions, the 
results of the work reported here are quite encouraging. First, the 
theoretical calculations predict the disappearance of the saddle shape, 
a phenomenon observed by many investigators. Second, the snap-through 
or appearance of two stable equilbrium states is predicted, another 
phenomenon observed by investigators. Finally, the transition from 
stable single-valued saddle solutions to stable cylindrical solutions is 
linked with a size effect. The figures show that both the thickness of 
the laminate and the length of the side determine whether the cylin- 
drical shape exists or whether the saddle shape exists. This investi- 
gator, as well as others, has felt a size effect exists in unsymmetrlc 
laminates and the model put forth here lends some credence to that 
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notion. 

While the predictions presented here exhibit all the Important 
features associated with unsymmetrlc laminates, several coniments are In 
order before closure. First, the solution presented here is a one-term 
Galorkln, or Rayleigh-Rltz, solution. Thus the solution Is, as with all 
one-term Galerkin solutions, over-constrained, This deficiency can be 
remedied by using more terras, and hence more generalized coordinates, in 
the assumed functional forms for u®, v®, and w. Modifying the current 
approach this way would probably change the numerical values associated 
with each solution branch of figs. 2 and 3 but not the main features 
of bifurcation and triple-solution. Second, the effects of moisture 
absorption, viscoelastic relaxation, or any other mechanism that alters 
the interital stress state of the laminate Is felt to be important for 
laminates V'jith lengths near the critical length. Alteration of the 
Internal stress state most likely Influences the numerical value of the 
critical length. Thus a laminate sized just above the critical length 
could, v^ith time, actually be sized just below the critical length due 
to any of the above mentioned time-dependent effects. Related to this 
is the fact that unsymraetric laminates sized near their critical length 
could exhibit "strange" behavior, requiring practically no force to snap 
them from one shape to another. More than likely, these multiple shapes 
would be some barely stable combination of shallow cylinders and shallow 
saddles. 

Finally, the analysis presented here is based on symmetric curing 
(symmetric about the z“0 plane) of the laminate and, as noted above, the 
lack of any effects to alter the internal stress state over and above 
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that due to temperature change during the cool-down from curing. The 
result is that for lengths greater than the critical length, two sira- 
illar shapes are possible. Each of these shapes has the same possi- 
bility of actually occurring, However, any external perturbation which 
la unsymmetrlc with respect to the midplane will cause one or the other 
of the two possible cylindrical shapes to be favored. Such a pertur- 
bation could be unsymmetrlc curing (due to unsymmetrlc cooling) , unsym- 
metric moisture absorbtlon, or any other process over which our control 
is limited. Thus in reality it may require, for example, more force to 
snap the cylinder of fig, Ic to the cylinder of fig. Id than it does to 
make the reverse snap. 
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C3 - A 3 ,lJ /1280 
C 4 - Ai2l2/48 
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The coefficients Aj^ 2 » ^22’ ®11 ^22 familiar defini- 

tions associated with classical lamination theory. 
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Fig, 1 Laminate Shapes: (a) at the elevated curing temperature, and 

at room-temperature, (b) a saddle shape, (c) a cylindrical 
shape, (d) another cylindrical shape. 
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Fig. 2 Room-uemperature shapes of square [02/902]„ T300/5208 graphite' 
epoxy laminates 
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Fig. 3 Room-Temperature shapes of square [0^/90^]^ T300/5208 graphite- 
epoxy laminates 
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